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Abstract. We show how to fix the renormalization scale for hard-scattering exclusive processes such as 
deeply virtual meson electroproduction by applying the BLM prescription to the imaginary part of the 
scattering amplitude and employing a fixed-t dispersion relation to obtain the scale-fixed real part. In this 
way we resolve the ambiguity in BLM renormalization scale-setting for complex scattering amplitudes. 
We illustrate this by computing the H generalized parton distribution at leading twist in an analytic 
quark-diquark model for the parton-proton scattering amplitude which can incorporate Regge exchange 
contributions characteristic of the deep inelastic structure functions. 

PACS. 11.55.Fv Dispersion relations ll.10.Gh Renormalization 


1 BLM renormalization scale setting 

A typical QCD amplitude for an exclusive process can be 
calculated as a power series in the strong coupling con¬ 
stant 

A = + AWa,(ti 2 ) + A^a 2 s (fx 2 ) + ■■■ (1) 

The renormalization scale /i of the running coupling in 
such processes can be set systematically in QCD without 
ambiguity at each order in perturbation theory using the 
Brodsky-Lepage-Mackenzie (BLM) method 

The BLM scale is derived order-by-order by incorpo¬ 
rating the non-conformal terms associated with the (3 func¬ 
tion into the argument of the running coupling. This can 
be done systematically using the skeleton expansion 
The scale determined by the BLM method is consis¬ 
tent with (a) the transitivity and other properties of the 
renormalization group (b) the renormalization group 
principle that relations between observables must be in¬ 
dependent of the choice of intermediate renormalization 
scheme j713| . and (c) the location of the analytic cut struc¬ 
ture of amplitudes at each flavor threshold. The noncon- 
forrnal terms involving the QCD f3 function are all ab¬ 
sorbed by the scale choice. The coefficients of the pertur¬ 
bative series remaining after BLM-scale-setting are thus 
the same as those of a conformally invariant theory with 
(3 = 0. In practice, one can often simply use the flavor 
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dependence of the series to tag the nonconformal (3 de¬ 
pendence in perturbation theory; i.e., the BLM procedure 
resums the terms involving n / associated with the running 
of the QCD coupling. 

Non-Abelian gauge theory based on SU(Nc) symme¬ 
try becomes an Abelian QED-like theory in the limit Nc —> 
0 while keeping a = Cfol s and ne = n e s/2CF fixed 0. 
Here Cf = {Nq — l)/2 Nc- The BLM scale reduces prop¬ 
erly to the standard QED scale in this analytic limit. 
For example, consider the vacuum polarization lepton- 
loop correction to e + e~ —> e + e _ in QED. The amplitude 
must be proportional to a(s) since this gives the correct 
cut of the forward amplitude at the lepton pair threshold 
s = 4to 2 . Thus the renormalization scale [Ir = s is exact 
and unambiguous in the conventional QED Goldberger- 
Low scheme m- If one chooses any other scale pL 2 R ^ s, 
the scale h 2 r = s will be restored when one sums all bub¬ 
ble graphs. The BLM procedure is thus consistent with the 
Abelian limit and the proper cut structure of amplitudes. 

2 Difficulties in using the mean value 
theorem to set the BLM Scale 

The BLM scale at leading order has a simple physical 
interpretation: it is identical to the photon virtuality in 
QED applications and the mean gluon virtuality in QCD 
when one uses physical schemes which generalize the QED 
scheme such as the pinch scheme |2in~n and the a v O 
E| scheme defined from the QCD static potential. The 
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number of flavors active in virtual corrections to a given 
process is evident from the BLM scale choice: the BLM 
method sets the renormalization scale so that flavor num¬ 
ber is changed properly in any renormalization scheme, 
including the MS scheme 12031 - 

In effect, the BLM prescription identifies the renor¬ 
malization scale and the gluon virtuality by eliminating 
the dependence on the number of flavors from the 0(a s ) 
(expanding a s itself) and 0(or s ) nonconformal terms in 
the perturbative amplitude. Typically, a QCD amplitude 
involves an integral over the momentum running through 
the gluon propagator. Therefore the argument of the cou¬ 
pling, if taken to be the momentum flowing through the 
gluon, varies through the integration phase space. A mean 
Q 2 can be extracted from the integral if the mean value 
theorem (MVT) of integral calculus can be applied. The 
essential requirement for the applicability of the MVT is 
that the function being evaluated at its mean value has to 
be continuous through the interval and the weight function 
has to be Riemann-integrable; this includes weight func¬ 
tions bounded and continuous in the range of integration. 
However, these necessary conditions are not a property of 
a principal value integrand associated with a pole which 
appears, for example, in Compton scattering and deeply 
virtual meson electroproduction (DVME); as we show ex¬ 
plicitly in Section IP below, the MVT does not apply for 
these amplitudes. 

It has been recently pointed out in Ref. m that the 
MVT prescription for determining the BLM procedure 
fails for amplitudes which are genuinely complex, that is, 
display non-vanishing real and imaginary parts. The au¬ 
thors then argue that there is no guarantee that the BLM 
prescription will yield the same answer for both parts of 
the amplitude. Worse, in the particular example that they 
study, p and tt\ electroproduction, the scale obtained from 
the real part becomes discontinuous (zero to infinity) at a 
particular kinematical point due to a divergence in some 
of the intermediate functions. 

In this paper we note that in a quantum field theory, 
the real and imaginary parts of a scattering amplitude 
are not independent but are constrained due to causality, 
locality, and Lorentz invariance. This manifests itself in 
the form of the dispersion relations traditionally used in 
meson photoproduction mm to link both parts of the 
amplitude. By examining a simple example in the next 
subsection, we show that the correct prescription for find¬ 
ing the BLM renormalization scales is to first fix the scale 
in the imaginary part, and then subsequently, the real 
part can be obtained by means of a dispersion relation. 
The BLM scales for the real and imaginary parts are thus 
not generally equal. We perform an explicit calculation 
for longitudinally polarized vector meson electroproduc¬ 
tion at nonzero skewness in Section I T2\ 


2.1 Example of the failure of the Mean Value Theorem 

To illustrate how the MVT can fail for an amplitudes 
which contains a pole, consider the following simple in¬ 


tegral: 

i[f](y) = t m dx M . (2) 

J -1 x-y + ie 

which is a functional of / and a function of y. The MVT 
would state that there would exist a certain x(y) in the 
interval (—1,1) such that 

i[f](y) = my))m(y) ■ (3) 


For / = constant the theorem holds trivially. But now con¬ 
sider a linear function f(x) = Ax + B , with A, B arbitrary 
real constants. Substituting Cauchy’s principal value for 
distributions (the imaginary part arising for y £ (—1,1)) 


1 

x — y + ie 


= PV 


1 

x-y 


inS(x — 


one easily finds 

I[f](y) = 2 A+(B + Ay) ^log 
= 2A + f(y)I[l](y) . 


y)> 



(4) 

(5) 


The MVT thus holds for the imaginary part, but it fails 
for the real part by the term 2A. Even in the instance 
where / could be chosen so that a certain x would satisfy 
the theorem for the real part, there would be no reason 
for x to be the same for the imaginary part. 

In this example, the function I[f]{y) is cut in the com¬ 
plex y plane in the interval [—1,1] of the real axis. The 
function satisfies an unsubtracted dispersion relation: 


nm = f 

2t n J c, 


dy'l[f){y') 


Rei[f]{y) = 


cut y 

PV 


W - 


( 6 ) 


y-y J 


dy'lml[f\(y'), 


and thus the discontinuity across the cut is sufficient to 
reconstruct the whole function. In this trivial example the 
discontinuity is 


I[f]{y + *e) - I[f]{y - it) = -27 rif(y) (7) 

and the mean value “scale” for the imaginary part is x(y) = 
y. Once this x has been chosen, the full amplitude is easily 
reconstructed from the dispersion relation. 

This procedure is not restricted to linear functions. For 
an arbitrary polynomial one can expand around y , 


/w = /(y)+E^T^(^-2/r 

n\ 


n =1 


( 8 ) 


we have now 

N r(n)( \ 

m ( y) = mm (y)+ E ((1 - y)" - 1 - y ) n ) ■ 

n= 1 

Again the second term is “unexpected” and a naive ap¬ 
plication of the mean value theorem to the real part fails. 
This additional term is a sum of binomials of increasing 
degree, which in the limit N —> oo, can be used to con¬ 
struct any entire function. Since it is analytic, it provides 
no contribution to the discontinuity across the cut, and 
the dispersion relation in Eq. © still holds. 
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2.2 Meson electroproduction 


We now turn to the problem of determining the renormal¬ 
ization scale for vector meson electroproduction 7 *p —> 
V°p', the critical example studied in Ref. (21 ■ We will 
employ the same kinematics, frame choice, and conven¬ 
tions as in Ref. El. It will be useful to note here that 
in the asymmetric frame employed, the skewness variable 
plays the role of Bjorken’s x, and 

C= Q*+ Q a -M% ] Q 2 = s + t + u ~ 2M N ■ ( 10 ) 

The skewness is in the interval 


Ce 


°’ 2M 2 



( 11 ) 


As a specific example of meson electroproduction we 
consider the production of a p° meson by a longitudinally 
polarized virtual photon 7 / with large virtuality Q 2 = 
-q 2 GOT- The relation of the differential cross section 
to the generalized parton distributions is well known m 

= _1_I V \M L \ 2 (12) 

dt 167r(s-M^i/2( Sr(? 2 )M 2 )2 LI 11 J 

where A 1 / 2 is the Kallen’s function, 

A 1 / 2 (a, 6 , c) = \J a 2 + b 2 + c 2 — 2 (ab + ac + be) (13) 


and the amplitude to leading order in a s and Q , ignoring 
the E Generalized Parton Distribution (GPD), is 

iMjji = -ie^ 1 47ra s (/r 2 )^/A(P)7 + ClA'(P , ) 

* L dZ<Pp ^4zP+ J c f xH pl( x ’^ ( 14 ) 

1 1 
X - 7 ---- _ 

x — Q — it x + le 

We use a model for the GPD H in this expression specified 
below in section^ Much useful information on GPD’s has 
been collected in E2 Notice that at large s one should 
also take into account double gluon exchange between the 
nucleon and the photon projectile (exiting the reaction as 
a meson). Substituting now 

Pa(P)7 + Pv (P') = 2 P + S X y (15) 

and employing the asymptotic value for the vector meson 
distribution amplitude m 


Once the BLM scale is fixed, as in Eq. iPl) below, the 
real part of the amplitude to leading order can be accessed 
by a dispersion relation 

1 f°° r 1 

Re Mp 0 = — Im Mh iPV —— 

n J(M N +m p ) 2 L S — 

A few remarks are in order: the ^ factor in Eq. G3 makes 
the dispersion integral in Eq. m convergent at fixed / 
(the Bjorken limit). This is the kinematical case tradition¬ 
ally considered in meson electroproduction m Further, 
the threshold for the integral is irrelevant if one consid¬ 
ers the large Q behavior of Mj' 0 . The use of a dispersion 
relation at fixed positive skewness £ induces a left cut in 
the complex s plane due to the photon becoming timelike 
there. This and the usual t-channel left cut are consis¬ 
tently neglected in eq. El as they are 1 /s ( 1 /Q 2 at fixed 
C) suppressed. At last and most important, the 1/Q factor 
in Eq. El relative to the corresponding Compton ampli¬ 
tude, makes the Compton cross section fall much slower in 
Q than its meson electroproduction counterpart. This is 
due to the extra gluon necessary to construct the LO dis¬ 
tribution amplitude of the emitted meson and is a striking 
manifestation of the J = 0 fixed pole which appears in the 
Compton amplitude due to the quasi-local coupling of two 
currents to the propagating quark. m 

Let us now turn to setting the BLM renormalization 
scale for meson electroproduction. We shall disregard the 
p evolution of H as an irrelevant complication for our ap¬ 
plication. Since the leading order is a Sl and the NLO am¬ 
plitude is available, we can follow m One identifies the 
NLO contributions from the vacuum polarization (which 
can be picked up from its dependence on Nf) and imposes 
the BLM scale-fixing condition on the imaginary part of 
the amplitude 



r 1 r 1 d x 

lm Jo dz ® p ^ J ( x 1 -Q/o H p°S X '^- T = {) ( 19 ) 


with the MS expression 


T = 


1 


z(x — ie) 

1 


- log (- log 


1 


V c 

5 


Q 2 


( 20 ) 


(1 — z) x — C + 

Q 2 


- log 


(l-z)(a;-C) 

c 


-log 




Following E] we first notice that the 2 integrals are trivial 
when employing the meson distribution amplitude in Eq. 

El . 


<P p (z) = 6z(l - z)f p 
with f p = 0.216 GeV we obtain 


(16) 


dz<P p (z) 


Im Mp 0 = 7) ~^ 47 r 2 ^ (VblmW aa ' ^ (17) 


dz<P p (z) 


logz = 9/p 

z ~ 2 

log(l - z) _ 9/p 

2 


1-2 


= ( 21 ) 


f dz@ p (z)- = 3fp 
Jo z 
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which allows us to find simpler expressions than those in 
m- Solving for p? yields the BLM scale 

p 2 = Q 2 e~ M/M (22) 

in the MS scheme where 


a 2 = 7r[J?(o,c,t)-tf(C,C,t)] 

'19 


Aj = 


6 


+ log C ) A 2 - 7 T (-H (C, C, t ) log(l - C) 


Jc-i x 

~ l dx 


(23) 


H( x,C,t) - H(C,C,t) 
x — C 


We also note that if Eq. GH> is Nf independent then Eq. 
m is also Nf independent, and the goal of the BLM 
scale fixing procedure is achieved simultaneously for the 
real and imaginary parts of the amplitude thanks to the 
dispersion relation. 

We thus have a totally consistent scale setting proce¬ 
dure at the level of the amplitude. Choosing to fix the 
BLM scale at the level of the total cross section is there¬ 
fore an unnecessary complication. The reader may observe 
that since both Ai and A 2 are linear in H , the scale in 
Eq. 1221) is independent of the absolute normalization of 
H. This would change if the E or other GPD’s were in¬ 
cluded in the analysis. The same observation applies to 
the normalization of the p meson distribution amplitude; 
i.e., p is independent of f p in our approximation. Thus the 
BLM renormalization scales are effectively independent of 
the flavor of the meson produced. 


3 Covariant quark-diquark model for the H 
generalized parton distribution 

To complete the evaluation of this specific example, we 
need a physically inspired model for H. We recall the an¬ 
alytic model for virtual Compton scattering introduced in 
Ref. [25] which is the gauge invariant leading twist exten¬ 
sion of the “covariant parton model” for structure func¬ 
tions given in Ref. m ■ This model can also incorporate 
not only Reggeon exchange in the t channel, but also the 
g 2 — and t- independent J = 0 fixed-pole contribution from 
the local coupling of the two photons to the struck quark. 
The J = an = 1/2 Regge contributions to the structure 
functions are evidenced by the ^fx. behavior of the nonsin¬ 
glet structure function F 2 p (i’,Q 2 ) — F%(x,Q 2 } at small x. 
Since this model has not been widely used in the past few 
years, we find worth readdressing it in the modern context 
as a contribution to the current discussion on GPD’s. 

The central ansatz of the covariant parton model is 
to construct a model of the quark-parton scattering am¬ 
plitude, a vertex with two quark legs of momentum fc, k! 



Fig. 1. (a) Quark-proton scattering matrix, (b) Quark- 

proton u-channel diquark exchange, (c) Quark-proton t- 
channel pomeron/Reggeon exchange. The need for modeling 
Regge trajectories is avoided by requiring t -C — 1 GeV 2 , where 
a(t) -> -1E2 and u aM oc A, which can be reabsorbed in the 
quark-diquark model. 


and spin indices i,i' and two proton legs with momentum 
p , p' and helicities A, A'. Since the parton legs are not on- 
shell, this amplitude is a function of four different Lorentz 
scalars, that can be chosen as the three Mandelstam in¬ 
variants s pp = (p + k ') , tpp = ( k' — k) 2 = ( p' —p) 2 , u pp = 
(p— k ) 2 and k 2 . The squared momentum of the returning 
parton can be expressed as k' 2 = s P p+tpp+u P p — 2Mf f — k 2 . 
We will denote this amplitude by 




pp 5 ^ppi Upp 


,k 2 


(24) 


In terms of the DVME scattering kinematical variables x, 
An, C, t, and the integration variable k, we have the fol¬ 
lowing expressions for the parton-proton scattering matrix 
Lorentz invariants: 

Spp = (1 + x - C )P + (k~ + ) - ( k J- ^ An) 2 

Upp = ( P + k~ - M^)(x - 1 ) - 

tpp = t (25) 

k 2 = xP + k~ - kl 

k 2 = t + M 2 r ( 1 - x) + P + k~{x - C) + {M 2 n + Ai) 
x (l + —^ - (kn - An ) 2 


(hereafter the subindex pp can be dropped without confu¬ 
sion). To connect with the modern formulation of deeply 
virtual Compton scattering in terms of GPD’s, let us note 
PT| that 


/\ 

H 0 x , (, t) 7 + - E{x, C, t ) 2M ^ 
<* p+ v-/ 2 {P',\'\V(0)' 1 +*{y)\P,\)\ y+= 


-e~ ¥y ~' 2 
4 fc 


U{P)> 


y+=0, y ±=0 


/l 

APsizp* - 1 + ) / AK. s i{x _ 0 P+ - *'+) 

(27t) 4 (5 (4) (p + k' - p' + k)Tx\'-i'i[s, t , u, k 2 } (26) 
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where the Fourier modes of the fermion field are 

$(k) = J d 4 ye ik y*(y)6W( yi _)6(y + ) . (27) 

This allows one to construct the leading twist handbag 
contributions to the GPD's given a model for the parton- 
proton scattering amplitude, by integrating over parton 
transverse and (—) momentum, and contracting the spins 
with a 7 + matrix in the parton Dirac space. 

The simplest model for the parton-proton amplitude 
is a tree level diagram based on perturbation theory with 
a proton-quark-diquark vertex (see Fig. Since we are 
only interested in leading twist effects in the electropro¬ 
duction amplitude, we will ignore the axial vector diquark 
and the spin dependence of the vertex. The proton spin in 
this model, with an s-wave proton-quark-diquark vertex, 
is carried by the struck quark, so that 

k 2 } =T[s,t,u,k 2 }6\y5ii / (28) 




x<0 


0<x<t 



Fig. 2. Integration over the k = k P + variable in the com¬ 
plex plane for the diquark model. From the figure, H(x < 0) = 


0. 


This can be viewed as a derivative with respect to A of 
Eq. Id 111 : since it can be pulled out of the k integration 
as a parametric derivative, it does not alter the analytical 
structure of the DVME amplitude nor the placement of 
poles of T in the k~ plane which we now study. 

Consider therefore the model defined by Eq. and 
l!321> . We employ the variable k = P + k so that an in¬ 
verse power of P + comes out of the k~ integral and yields 


for the H GPD part, that yields an oversimplified yet prac¬ 
tical 

2 p+H(x,Ct) = ^ J -^^S{P + x - k + )T[s,t,u,k 2 } . 

(29) 

The amplitude in covariant perturbation theory is 


iT(s,t,u,k 2 ,k' 2 ) = ( ig(k 2 )) _ l _ ^ + .J ,ig{k 2 )) 

(30) 

Here A plays the role of a scalar diquark mass and g is the 
amplitude for the proton-quark-diquark coupling. The sta¬ 
bility condition dictates that the sum of the diquark and 
parton masses has to be larger than the proton mass to 
prevent the decay to a free quark and diquark pair. This 
feature is not apparent from Dyson-Schwinger models of 
the proton-quark-diquark vertex EBl, EH!, mm, where the 
diquark mass is of order 800 MeV and decreases simulta¬ 
neously with the quark mass at large momentum, but the 
Euclidean space formulation does not accommodate the 
decay. Except for this point, the work of these authors 
can perfectly be reused here as model input. 

The vertex functions g{k 2 ) are known from these fun¬ 
damental studies to fall with k 2 . One could guess 


g(k 2 ) 


9 

( k 2 — A 2 ) + ie 


(31) 


for A some constituent or running quark mass (that we 
will take fixed here at 0.5 GeV, with g = 55 and A = 
1.6 GeV). This is simply suggested by the observation 
that T is not a truncated amplitude at the quark legs, and 
therefore it includes their propagators. But since all the 
recoil t dependence in the form factor calculation needs to 
be present in the vertex g{k 2 ) and we know by the QCD 
counting rules that at large momentum transfer Fi(t) oc 
V 2 , the simplest workable model is 


g(k 2 ) 


_ g _ 

((fc 2 - A 2 )+ief ‘ 


(32) 


H( X ,c,t)= r r r dK-T[sxvk 2 \. 

Jo \^ 7r ) Jo Jo 

( 33 ) 

in place of Eq. nu and with k + , k’ fixed as above. The 
parton-proton amplitude is at this order an holomorphic 
function of k~ and has three poles, whose positions are 
depicted in Fig. [3 The diquark propagator yields a simple 
pole denoted as (1) and given by the condition it—A 2 Vie = 
0. The vertex functions yield two double poles respectively, 
denoted by (2) for k 2 —A 2 Pie = 0 and (3) for k 2 —A 2 +ie = 
0. These poles are respectively located at 


A 2 + - Mjy(l - x) 

«i = -i- 

A 2 + + ie 

Kr, = - 

x 

A 2 + c 2 + ie 
k 3 = -7- 

x — C 

where we have grouped 


(34) 


t + 1 — x) + (Mfj + A^_) 


x 1 + 


1-C 


- (kj_ - Aj 


(35) 


The k~ integral can now be immediately performed 
picking up the residue of the simple or double pole as in 
Fig. 01 We find 


H(xe(0,OX,t) = g 2 A s L^ 
*1 Hk’i-Aif) 


nOO 

/ |k x |d|k x | (36) 

Jo 

( x ~ 1 + 2 ( a: ~CA 
\ it - A 2 k' 2 - A 2 J 


evaluated at the pole k~ = k^ ■ The second term, be¬ 
ing proportional to x — £, does not give a contribution to 
the imaginary part of the electroproduction amplitude at 
leading Q 2 . We also have 
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Fig. 3. Two light-cone time-ordered diagrams for the hand¬ 
bag diagram with a covariant parton-proton scattering vertex 
based on a parton-diquark-proton vertex and second order per¬ 
turbation theory. Notice that the second diagram, present in a 
covariant model, gives rise to a non-vanishing H(x £ (0, (),(, t ) 
and to a well-defined imaginary part, by providing a finite part 
from the left needed for continuity at x = 

H GPD in diquark model 


Fixed ;=0.12 



Fig. 4. The t dependence of the Generalized Par- 

ton Distribution in the perturbative diquark model at fixed 
£ = 0.12. Since this model does not have a sea wavefunction 
component, it does not give access to the region £ — 1 < x < 0. 
However it is sufficient for the purpose of illustrating BLM 
scale fixing via a dispersion relation. 


H(xe ((,l),(,t)=g 2 A 8 



1 ~ C/2 

(2tt) 3 

1 


|kj_|d|k 

1 


' 2 - A 2 ) 2 (k 2 - A 2 ) 2 


(37) 


which gives a positive definite H function. 


4 Numerical results 

We now provide a numerical computation of the above 
model in order to illustrate the approach. The t depen¬ 
dence of the H GPD in this quark-diquark model is plot¬ 
ted in Fig. Q]. As can be seen, H falls rapidly with t at 
fixed The region £ — 1 < x < 0 cannot be accessed with 
this “valence” model. However, the fact that the vertex 


H GPD in diquark model 

Fixed t=-2GeV 2 



Fig. 5. £ dependence of the H(x,£,t) Generalized Parton 

Distribution in the perturbative diquark model at fixed t = 
— 2 GeV 2 . Observe that by performing different pole integrals 
for x > £ and x < £, obtaining continuity at x = £ is a non¬ 
trivial test of the model and computations. 


employed is covariant gives a feature absent in the con¬ 
stituent valence light-front quark model, that is, a non 
vanishing H(x = £, £, i) and continuity at x = £. This 
allows for an imaginary part in the electroproduction am¬ 
plitude in this minimal model. In Fig. [3 we plot the £ 
dependence of H. The sharpness of the peak in this graph 
is controlled by A, the diquark mass. The diquark model 
satisfies the simplifying property 

H(x < 0,C,i) = 0. (38) 


We then can reproduce the result of Belitsky and Muller 
for the BLM scale for the imaginary part of the amplitude; 
in our asymmetric frame Eq. |23|) their result is 


M 

a 2 


19 , 

= ¥ + !og 




dx H «' M , 

x-( 


(39) 


In Fig. E we display the resulting BLM scale as ex¬ 
tracted from the imaginary part of the amplitude. It can 
be seen that fi 2 BLM <C Q 2 , indicating that meson electro¬ 
production at practically accessible kinematics measures 
the coupling constant a s in the infrared regime where the¬ 
oretical considerations m lead us to expect a fixed point. 
For numerical purposes one can use 


47T 

] = ( 9 l°g((M 2 + Mq)/A 2 )) (40) 

with A ~ 0.2 — 0.21 GeV and Mq ~ 1 — 1.1 GeV at the 
BLM scale in the MS scheme. With this, eqs. mm can 
be immediately evaluated. This legacy form of the running 
coupling has been successfully applied to the prediction 
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BLM scale in Meson Electroproduction 

Fixed by imaginary part 


T 1 1 1 T 1 

1 T 1 1 1 


- 

• Total 

— contribution of integral over H 

- 



- 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

_1_1_1_1_1_1_1_1_ 

; 

_1_1_1_i_1_i_ 


p° L electroproduction 

t=-2GeV 2 0.3 



Fig. 6. BLM scale in p° L meson electroproduction at Q 2 = 
4 GeV 2 and t = —2 GeV 2 . The solid line represents the con¬ 
tribution of the integral in Eq. 


Fig. 7. As an illustration we present the model computation 
for p° L meson electroproduction at t = —2 GeV 2 and ( = 0.3. 
The fall-off with Q 6 at fixed £ follows from Eq. 1141 and 1121 . 


of electromagnetic form factors and other quantities (see 
overview in E21). Since the BLM scale that we determine 
is small, our calculations sample mostly a s (0). The value 
of this fundamental quantity, controlled in eq. m by the 
parameter Mo, is uncertain. With our parameterization, 
a s (0) ~ 0.9. 

As a final illustration we give the model’s prediction 
for the differential cross section in Eq. O- The result 
of the computation is plotted in Fig. [3 Since the per¬ 
turbative diquark model is only viable at sizeable t, we 
fix t = —2 GeV 2 . Consequently the cross section is much 
smaller than reported by other authors at t = ^minimum- 
The fall-off as a power law at fixed skewness is built-in 
into our formalism, and can be traced to Eq. d , <tT21 1 . 
(Notice that the DVCS cross section is enhanced at high 

o 2 

momentum transfer by a factor ^ 2 - and eventually crosses 

over and dominates over p° production in spite of the ex¬ 
tra a.EM suppression in Compton scattering [223 )• For this 
example we have chosen kinematics such that Q 2 /t is 2 
or more as a compromise between theory and experimen¬ 
tal uncertainties. From the theoretical point of view one 
would like to control the factorization corrections increas¬ 
ing this ratio. However a much larger Q 2 would be difficult 
for an experiment at an upgraded Jefferson lab or, due to 
statistical limitations at sizeable t, at larger s machines. 
Note that the largest uncertainty in our work likely comes 
from modeling the GPD H. The DVME Q~ 6 fall-off is a 
genuine QCD prediction m and consistent with the ex¬ 
perimental result of HERMES [32| at s = 5.4 GeV that 
yields an exponent —5.66(36), within one standard devia¬ 
tion. One should notice the subtlety that the Hermes data 
is quoted at fixed center of mass energy squared W 2 = s 
for the photon-nucleon system, whereas the predicted fall- 
off is at fixed scaling variable, here the skewness £. But 
since reference 1321 also displays a very weak W depen¬ 
dence, one can slide s in the Hermes data and use our 


eq. dU to obtain the same exponent at fixed skewness. 
The cross sections reported by this experiment, which are 
dominated by low t Regge processes, cannot be directly 
compared with our result. We look forward to future mea¬ 
surements with a larger —t reach. 

5 Conclusion 

We have shown that there is no a priori difficulty to set 
the renormalization scales which appear in hard hadronic 
scattering amplitudes that are genuinely complex. The in¬ 
dividual BLM scales for the real and imaginary parts are 
unambiguous and distinct, and they maintain the analyt- 
icity and dispersion relations of the full scattering ampli¬ 
tude. Just as in Abelian theory, the BLM method in QCD 
fixes the arguments of the running coupling to maintain 
the correct cut structure of the quark loop vacuum po¬ 
larization contribution at each flavor threshold. All the 
non-conformal contributions to the amplitude which are 
associated with the running of the gauge coupling are used 
to set the renormalization scales. The remaining perturba¬ 
tion series is then identical with a conformal theory with 
zero /? function. 

We have determined the BLM renormalization scales 
for the imaginary and real parts of the <y*p —> p' meson 
electroproduction amplitude, a process of current theoret¬ 
ical and experimental interest. We applied the BLM pre¬ 
scription to the imaginary part and employ a fixed-t dis¬ 
persion relation to calculate the scale-fixed real part. We 
have exploited the connection between generalized par- 
ton distributions and the parton-proton scattering ampli¬ 
tude in order to obtain an analytic representation of the 
helicity-conserving amplitude H in terms of a simple, yet 
physically inspired, quark-diquark model of the proton. 

It has been conventional to characterize the precision 
of pQCD predictions by using an arbitrary renormaliza- 
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tion scale in the MS scheme, such as /i^ = Q 2 , and then 
varying the scale over an arbitrary range, e.g., Q 2 /2 < 
fj? R < 2 Q 2 as a means to estimate the convergence of the 
perturbative series. However, the variation of the renor¬ 
malization scale can only be relevant for the nonconformal 
contributions to an observable, not the complete series. In 
fact, as we have stressed, there is actually no renormaliza¬ 
tion scale ambiguity since one must set the argument of 
the running coupling such that amplitudes have the cor¬ 
rect analytic cut structure at each quark threshold. This 
is done correctly in any renormalization scheme by using 
the BLM method. 
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